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SUI>1MARY 


A mathematical model of the Sikorsky C‘H-3G helicopter based on classical non- 
linear, "quasi-steady" rotor theory has been developed .it NASA Ames Research Center. 
The model has been validated statically and dynamically by comparison with Navy 
flight-test data. The model incorporates ad hoc revisions which address the ideal 
assumptions of classical rotor theory and improve the static trim characteristics to 
provide a more realistic simulation, while retaining the simplicity of the classical 

model . 


INTRODUCTION 


The Guidance and Navigation Branch at Ames Research Center is conducting research 
to improve helicopter IFR operations at remote sites and at high-density traffic 
areas. Much of the research is accomplished using a Sikorsky SH-3G helicopter (see 
fig. 1) to evaluate advanced guidance-and navigation concepts. Prior to flight test, 
new* concepts are developed on an off-line simulation or using a real-time piloted 
simulation. A requirement, then, exists to develop and validate an off-line math 
model of the SH-3G which can be adapted for real-time simulation. 

In recent years, NASA has developed several simulations of Sikorsky aircraft. 

In 1979, J. D. Shaughnessy of Langley Research Center developed a math model of the 
Sikorsky CH-54 helicopter for sling-load research (ref. 1). The rotor models used 
were based largely on an NACA report by F. J. Bailey (ref. 2), who related rotor per 
formance to only three varying parameters: the inflow ratio, the tip-speed ratio, and 

the rotor pitch. However, Bailey assumed uniform downwash which leads to underesti- 
mating the induced power by approximately 11% in hover, and 17% in high forward flight 
(ref. 3, p. 140). In 1980, Sturgeon and Phillips (NASA Ames) modified Shaughnessy s 
model to simulate the Sikorsky CH-53 (ref. 4). 

This paper documents a mathematical model of the SH-3G helicopter which was 
developed by modifying and adding to the existing CH-53 helicopter math model at 
Ames, and validated by matching flight data. The present model differs from the 
CH-54 and CH-53 models in that actual static performance as measured in flight test 
(ref. 5) is more closely matched by addressing the assumptions of uniform downwash, 
two-dimensional, blade lift curve slope, and fuselage flat-plate area as measured in 
the wind tunnel. In addition, the fuselage aerodynamics, equations of motion, and 
engine model are simplified. An improved trimming algorithm has also been implemented. 

Like the CH-53 math model, the SH-3G math model calculates nonlinear rotor aero- 
dynamics based on the "quasi-steady" assumption, i.e., there are no unsteady aerody- 
namic effects between time steps. The fuselage aerodynamics have been line.arized as 
much as possible and the engine and associated governor are modeled by a simple trans- 
fer function between the main-rotor rpm deviation from nominal and tne engine torque. 





Figure 1.- NASA SH-3G research helicopter. 


MATHEMATICAL MODEL 


The SH-3G math model consists of eight submodels: the equations of motion, the 

atmospheric model, the wind, the fuselage aerodynamics, the main rotor, the tail 
rotor, the engine, and the control system and rigging. A description of the 
coordinate systems used throughout, a general description of all the submodels, and 
a detailed description of each submodel follow. 


Coordinate Systems 


1. Earth axes, subscript e: Origin fixed on the earth's surface, x axis 

pointing north, y axis pointing east, and z axis pointing down into the earth. 
This coordinate system rotates with the earth with the z axis always po n ng 
toward the earth's center. 


2 . 

X axis 
dicular 
z axis 
fig. 2) 

3. 

X axis 
dicular 
fig. 3) 


Path axes, subscript p: Origin at the center of gravity of the helicopter, 

pointing along the earth relative velocfty vector, y axis pointing perpen- 
to the right of the earth relative velocity vector and parallel to the ground, 
pointing down and perpendicular to the earth relative velocity vector (see 

Body axes, subscript b: Origin at the center of gravity of the helicopter, 

pointing out the nose of the helicopter, y axis pointing to the right perpen- 
to the plane of symmetry, and z axis down in the plane of symmetry (see 
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4. Shaft axes, subscript s: 
Origin at the rotor hub, x axis 
rotated about the y body axis through 
the longitudinal shaft tilt angle, 6g, 
y axis rotated about the new shaft x 
axis through the lateral shaft tilt 
angle, (^g, and z axis down and paral- 
lel to the shaft (see fig. 4). 

5. Control axes, subscript c: 
Origin at the rotor hub, x axis point- 
ing toward the relative wind parallel 
to the swashplate, y axis pointing to 
the right parallel to the swashplate 
(perpendicular to the relative wind) , 
and z axis down and perpendicular to 
the swashplate (see fig. 5). 

6. Wind axes, subscript w: 

Origin at the center of gravity of the 
helicopter, x axis pointing into the 
relative wind, y axis rotated about 
the z axis by the sideslip angle, 6, 
and z axis rotated about the y 
wind axis by the angle of attack, a 
(see fig. 6). 


General Model Description 

The SH-3G helicopter simulation 
contains the following submodels: 

1. Equations of motion: This 

submodel calculates the position, 
velocity, acceleration, attitude, 
angular velocity, and angular accelera- 
tion from the forces and moments pro- 
vided from other submodels. 

2. Atmospheric model: Atmosphere 

pressure, temperature, density, and 
dynamic pressure are calculated from 
the 1962 standard atmosphere. 


N 



Figure 2.- Path axes. 



3. Wind model: Turbulence and Figure 3.- Body axes, 

steady-wind components are generated 

in this submodel. The random turbulence conforms to the Dryden spectral model. 

4. Fuselage aerodynamics model: The fuselage aerodynamics model determines the 

lift, drag, and side forces, as well as the pitching, rolling, and yawing moments as 
functions of the fuselage angle of attack, the fuselage sideslip angle, and the 
dynamic pressure. In the interest of simplicity, the functions are all linear or 
polynomial functions of sinusoids. 
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Figure 7.- SH-3G math model block diagram* 

Equations of Motion 

The SH-3G equations of motion are a simplified version of SMART, a standard sub- 
routine for simulation at Ames and documented i, reference 6* SMART converts forces 
and moments from body axes to earth axes, integrates in earth coordinates to deter- 
mine the earth relative velocity and position and converts the velocities back to body 
axes . 

The SM-3(^ version of SM/\RT deletes the small earth Coriolis effects and uses 
equations l:or the standard atmosphere instead of a table. A more convenient initial- 
ization has also been implemented (see appendix A). 

Ihe body-axis torces generated in the fuselage aerodynamics model, the main-rotor 
model, and the tail-rotor model are summed to produce the total body-axis forces 
acting on the lielicopter. 

the total body-axis forces are related to tlie earth— axis forces by the familiar 
Euler ancle tuitions: 
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The trigonometric functions "sine" and "cosine" will be abbreviated "S" and "C 
throughout this paper. 

Solving for the earth relative forces: 
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If we neglect the earth Coriolis accelerations, but keep the earth centripetal 
accelerations, we find the earth’s relative accelerations to be: 

AA / \ Ae'< C(LAT)S(WT)\ 


F + F 
D C 


R C^(LAT) 
e e 


where Fr is the gravitational force equal to the weight of the aircraft, m is the 
mass of the aircraft, Rg is the radius of the earth, and is the rate of angular 

rotation of the earth. The aircraft velocity magnitude relative to the earth has 
been neglected compared to the inertial velocity due to the earth s rotation. 

The earth relative velocities are found by integrating the earth relative accel- 
erations. A second-order Adams-Bashford predictor algorithm for integrating is used 
(as in ref. 6) yielding: 


+ 3x V, 


where DT is the integration time step, n is the present value, (n-1) is the pre- 
vious value, and (n-fl) is the next value. 
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Earth location is now found by numerical integration using a modified 
method: - t 


LAT I + I LAT I + I hAl' 


where ALT is the altitude above sea level. 

1 • 4^1 ^r> rVe lir in earth coordinates can be calculated from 

The velocity with respect to tt.e air in eaiLn 

the vector equation: 
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where Va = ^0^+ ^ "•''^g^^tude of the air velocity and \%\ 0. 

If the x-z plane is a plane of symmetry, the body-axis, angular accelerations 
(from ref. 6) are; 


(C„P + C„R)Q 


' 8 ^ ' '"Tl 


'p\ (CiR + C;,P)Q C 3 0 C^ /O 
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where L, M, and N are the total body axis moments generated by 

the fuselage aerodynamics and . . • C^q are inertial coe icien s as 


Co = (Wzz - ^xz^"' 


Ct = C„[(Iyy - - 1^,1 

^2 ^ Cq^xz^^xx " ^yy ^ ^zz^ 
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3 0 2Z 
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0 XZ 
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6 7 xz 


c, = r 


Cp = C()[(I^J^ - Iyy)^XX ^XZ^ 

Cg “ Co^xz^^yy ” ^zz ~ ^xx^ 

Cl' = Colj^x 

where Ixx* St’ inertia. 

Integrating the body angular accelerations using the second-order Adams-Bashford 
uredictor algorithm yields the body-axis angular velocities: 



oWGiNR'- Pf;®, 
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The body axis acceleration vector is (see ref. 6): 


( i .:0 



(14) 


Neglecting the earth* s angular velocity (7.2722*E-5 rad/sec), the Euler angular 
accelerations are: 



f(Q S(f) + R C(J))/C0^ 
Q C<{) - R S(() 

P + SG 


(15) 


Integrating the Euler angular accelerations using tha modified Euler r!]gurithm: 




(16) 


Atmospheric Model 

The pressure, temperature, and density of the atmosphr q have important effects 
on the aerodynamics of any aircraft. 


Since helicopters operate well below the speed of sound, conditions can be 
assumed subsonic. Most helicopters also operate in the troposphere (below 36,089 ft 
in the 1962 standard atmosphere) so for this case the atmosphere can be modeled by 
just two equations: 


T. = 1 - 6.875 E-6 x ALT 


5.256 


P. = T. 


( 17 ) 


where T.^^ and are the temperature and pressure ratios relative to standard 

sea-level vahies and AT/I is the altitude above see level in feet. 

For a diatomic gas such as air, the total temperature and pressure ratios are 
given by: 
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T,- « ] + 0.2 M' 


P ;n: *1* ^ ^ 

J .1- y. 


/here M is the Mach number. 

The. ambient conditions are now found as: 

T = AT + To * T . 

'^r Htd 


P “ I’ar * '’std 


p = P/(T * 3088.8) (P in psf, T in Kelvin) 


V . = 65.76 
sound 


(T in Kelvin) 


where Tgtd ^^id Pstd are the standard sea-level values of pressure (21 16 ,2 
temperature (288.16 K) • 

The total pressure and temperature can now be found from equation (18) 

\ot = "r - I 

P = P A P I 
tot r / 

The compressible and incompressible dynamic pressures are: 

q = 0.5 pV^ \ 


psf) and 


q = P - P 
tot 


Finally, the equivalent and calibrated airspeeds: 


= '^'f^std * 

= l/5{[l + 


cal sound 


Wind Model 

The wind model JLncli'des steady winds and turbulence conforming to the Dryden 
spectral model. The" wind is assumed to be a frozen field of turbulence drifting at a 
mean wind speed relative to the earth. The turbulence is defined by a characteristic 
high-frequency cutoff wavelength and a standard deviation vhich is a function of 
altitude and intensity category, i.e., light, moderate. The cutoff wavelength is 
determined by the rotor diameter, as the rotor is assumed to be completely engulfed 
by each change in the wind from time step to time step. With this assumption the 
turbulence is generated by simply adding a random component to the velocity vector in 
body axes, rather than an elaborate integration across the rotor disc. 






The changes in the velocity vector enter the force- and moment-generation svib- 
routines and affect the accelerations of the lielicopter model. Angular accelerations 
are simulated by simply adding random increments (conforming to the Dryden model) to 
the angular velocity components. 

All random turbulence effects are removed during trim; however, the mean wind is 
included. The random velocity changes as a result of turbulence would make it 
impossible to trim. 

The wind- model used in the SH-3G simulation is identical to that of SMART (see 
ref. 6). A background description can also be found in reference 7. 


Fuselage Aerodynamics 

The fuselage aerodynamic forces and moments are presented as functions of the 
angle of attack, sideslip angle, and dynamic pressure. The forces and moments are 
first presented in coefficient form as functions of the two aerodynamic angles. The 
coefficients are then multiplied by the dynamic pressure and transformed to body axes. 

Data for the fuselage aerodynamics model are taken from two sources: a trainer 

math model of the SH-3H (ref. 8) and a Sikorsky Engineering Report (ref. 9), both 
prepared under Navy contract. Reference 8 contains equations for the various forces 
and moments and reference 9 contains actual wind-tunnel test data. 

The fuselage aerodynamics submodel has been greatly simplified from the CH-53 
subroutine on which it is based. The philosophy has been to use linear or trigonomet- 
ric functions to approximate the wind-tunnel data when simple equations are not 
already available. Model accuracy is maintained for small values of the angle of 
attack and sideslip angle, but no attempt is made to fit the data exactly for large 
angles. The large angles are, generally, only possible at low airspeeds which means 
small fuselage forces and moments compared to the rotor aerodynamics. Further 
sophistication of the model for large aerodynamic angles is therefore not justified. 

Fuselage damping moments have also been neglected. This does not substantially 
harm the fidelity of the simulation because the simulation will almost always be used 
with the automatic control system on to compensate for the lack of natural damping. 

The first parameter to calculate is the effective angle of attack of the fuselage 
accounting for the main rotor downwash. This local angle of attack is only used in 
finding the fuselage aerodynamics. The rotor downwash factor (from the CH-53 model) 
is; 


(23) 

III 

where Ct is the main-rotor thrust coefficient, is trie main rotor inflow ratio, 

in) m 

and Un) is the main rotor tip-speed ratio. 

The fuselage local angle of attack is: 


u 






( 24 ) 


wliere 


^kf 


is taken as 0.5 from reference 9. 


1 1 








S?Toor 

The following force and moment coefficients, in fuselage wind axes, have the 
dimensions of square feet. The aerodynamic angles are in degrees. 

The side-force coefficient is an approximation to wind-tunnel test data from 
reference 9: 


% = -73 


1 3 i 50 


= -350 sign(3) |3l > 50 ) 

The lift-force coefficient is a strong function of both the angle of attack and 
the sideslip angle. The angle of attack component is taken from the Navy trainer 

model (ref. 8): 

C = 10 + 410 S(a^) (2( 

The lift-force coefficient, as a result of sideslip, is a linear approximation 
of wind-tunnel data from reference 9: 


= -4 3 


= -40 


|3l < 10 
10 < |3| <40 


= -7.6 X 1 3| + 344 40 ^ I 3| < 65 


= -150 


3l > 65 


The drag-force coefficient is also a strong function of both angle of attack and 
sideslip angle. The angle-of-attack component is again taken from the Navy trainer 

model (ref. 8): 

c = C + 324 S^(ccj^ + 2) (28 

^Di 

The flat plate area, was determined ad hoc by matching the torque at 

90 knots. The value finally used was 44 square feet. 

The drag force coefficient because of sideslip is a sinusoidal approximation to 
Sikorsky wind-tunnel data from reference 9 * 

= 500 s'CS) (29 


Jhe following moment coefficients are with respect to body axes and have dimen- 
sions of cubic feet. They are all approximations to Sikorsky wind-tunnel test data 
from reference 9, except for the pitching moment due to angle of attack which is from 
the Navy trainer model. 

The rolling moment due *:o sideslip is given by: 


C{,j. = -4.5 3 


|b1 < 20 


-90 sign(3) ! R I >20 


Tho pitching'-momcnit coefficiont haw two compommts: ono from angle of attack 

and another from .sideslip. The component due to angle of attack is given by: 

170 - 1950 S((Xp) 
ml 

The component due to sideslip is 

= -175 S(4 ^31) 
m? 

The yawing moment due to sideslip is given by 

= -400 S(4.5 3) |3l < 40 

n 

- 80[3 - 40 sign(3)] 40 < |3l < 90 
= 4000 * sign(3) |3| >90 

The fuselage forces are found by multiplying the force coefficients by the 
dynamic pressure which is given by 

q = 0.5 pV^ 



The forces are then given by: 



In wind axes: 



The wind-axis forces are transformed to body axes a.s follows: 
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The fuse] ago. moments, in body axes, arc found by multiplying the dynamic pressure 
times the moment coefficients. 


T - q 
a M , ^ 


m \ mi m 2 / 


1 


(38) 


T = C q 
a m ^ 
n n 


Main-Rotor Model 


The main-rotor model is adapted from the CH-53 version (see ref. 4), which is a 
classical nonlinear Bailey model assuming uniform inflow. An empirical method to 
account for nonuniform inflow using flight-tesf data has been included in the SH-3G 
rotor models. 

The main-rotor and tail-rotor models are based on the following assumptions: 

1. Compressibility and stall effects can be neglected. 

2. Lag effects can be neglected. 

3. Only the first harmonic motion of the rotor blades is important. 

4. Blade coning and flapping angles are quasi-static. 

5. Any wind or turbulence emerses the entire rotor disk at once. 

The following discussion is paraphrased from reference 4 (for completeness) with 
the exception of the torque equation. 

The airspeed of the entire rotor disc is assumed to be that of the rotor hub. 

The airspeed at the rotor hub is calculated in shaft axes using the helicopter 



where is the shaft pitch angle; <))s is the shaft roll angle; xj^, yj^, are the 

hub coordinates in body axes; and U^, V^,, and Wj, are the body-axis velocities. 

The airspeed at the hub is transferred into control axes using the rotor orien- 
tation angle: 


S = arctani 
r 


/V + Ai W \ 
I s ^s s 1 

\U + Bi W y 
Vs s' 


( 40 ) 


which is obtained using the definition of control axes: = 0. 

Using small am'le approximations for the main rotor cyclic control inputs 
(swashplate angles), A,^ and Hig* 
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The Bailey method requires three parameters on encli iteration: tip-speed ratio, 

inflow ratio, and the rotor pitcli. The tip-speed and Inflow ratios are defined as: 


1. = lV/i?R 

\ = Wp/t'R - V 


(A2) 


where i: is the angular velocity of the rotor, R is the blade radius, and v is 

the induced inflow ratio. 


The induced inflow ratio, v, is found by filtering the steady-state value, an 
idea first used by Schaughnessy (ref. 1). The resulting first-order, nonlinear 
differential equation is more stable numerically than algebraic calculations. 


'v 


') 


(A3) 


The thrust coefficient is C-j^' and is an empirical time constant to simulate 

the lag associated with inflow changes. 

The rotor thrust coefficient, C.-p, and the coning angle ay (see fig. 8) are cal- 
culated by equations from ref. A. Tbese equations are simplified versions of those 
found in Bailey's NAi'A report (ref. 2). 


(AA) 


where e is the solidity ratio, ^ is a blade mass constant, B is the tip loss 
factor, t'y is the rotor pitch, I'j is the rotor twist, and a is the blade-lift 
curve slope. The blade-lift curve slope is usually taken as 3.73/rad based on two- 
dimensional, wind-tunnel data and the resulting thrust coefficient is usually opti- 
mistic. for the present simulation, the blade-!ift curve slope was determined 
empirieallv by decreasing its value until the simulation collective stick position in 
hover, matclusl flight -test data from reference 3. The final value was 5.2/rad. 

I'he calculation of the rotor flapping angles (sei> fig. h) r*'quires the fuselage 
angtilar vi'locity it control axes. This requires the following two t ransformat ions: 
from bodv to shaft axes and from shaft to control axes. 





AXIS OF NO AXIS OF NO 



Figure 8.- Coning angle. 




The flapping angles a^ and are calculated 

in control axes by formulas from reference 1: 



where 75 is the blade pitch at 75 % of the rotor radius, i.e., 




G 


I 7 b 


n 

0 


+ 0.75 0^ 


The rotor drag force in wind axes is given byt 

H = T*a' 


where a' represents a lift-induced, tilt back angle of the rotor thrust vector and 
Is given by (ref . 1 ) : 


1 


(47) 
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The Qj, term accounts for fuselage pitch rate. 

The torque equation is taken from page 151, equation (5.27a) in reference 3: 


C = ^ (1 + K X 3p^) + C„(K ,u - \ - ya') 
q 8 sw ' T nu 


(48) 


where is a constant to approximately account for spanwise flow in forward flight, 

similarly accounts for nonuniform inflow, and 6 is the blade profile drag 
coefficient and K^u were determined as follows: 

= 1.57 (attributed to Stepniewsky in ref. 3, p. 151) 

K =0.15 (empirically determined by matching the simulation torque in 
hover to flight-test data from ref. 5) 


The Bailey report (ref. 2) gives 6 as: 

6 = 6() + 


(49) 


where Uj. = 6Cx/(oa) which is the average blade angle of attack in hover. For the 
NACA 0012 airfoil and calculated from reference 10 = 0.0078, 6 ^ = -0.0090, and 
6.. = 0.2987. 

The torque is calculated from the definition of torque coefficient as: 


Q = bcR’p(Rt')^ — ^ 


(50) 


where b is the number of blades and c is the blade chord. 

The rotor side force is defined in terms of the rotor-side-force coefficient: 

C 


J = bcR, (R;') 


0 


(51) 


where 




I apii\ + a^biU - aoa,Vi^ + ^ a^^a^ 


m 

- (I ya„ - j b, - \ y 


Tills last equation is from reference 4 and was used on the CH-53. 


(52) 


I'he rotor forces are transformed from control axes to body axes by using the 
transposes of the matrices in equations (45). First transform control to shaft axes: 


17 


Of POOR 




Then transform from shaft axes to body axes: 




( 53 ) 


(54) 


The rotor cyclic controls enter the main rotor model through the cyclic flapping 
angles. From reference 4, these angles in shaft axes are given by: 



(55) 


Pure pitching and rolling moments about the hub are generated by the flapping 
hinge of^fts! A^ure yawing'moment is generated by the engine torque. In shaft axes 

these pure moments are given by (ref. 4): 



(56) 


where e Is the flapping hinge offset, is the blade mass moment , and Qe la 

the engine torque. 

The total moments of the main rotor in body axes result from the 
(converted to body axes) plus the rotor forces acting through the coordinates 
rotor hub relative to the center of gravity: 



CO , 


s S ■ S a S 

0 C<J) -S({> 


-S*! CO CO C<) 

S ,S S S S 



0 


-z, 


^h ° 


~yx 


^h 

0 



(57) 


where xh, Yh. and are the coordinates of the rotor hub relative to the center of 

gravity. 

The maln-rot.ir angular acceleration ia found by summing the torques at the hub: 


= - ^a - ^t • ^^/'m 


(58) 
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where Q^- is the tail -rotor torque, is the ^ear ratio between the main rotor and 

tail rotor (G^ > 1), and is the moment of inertia of the rotor blades and hub. 


Tail Rotor Model 

For the most part, the tail-rotor model calculates forces and moments in the same 
manner as the main-rotor model. The tail rotor does not have cyclic-control angle 
commands. The tail rotor does, however, have a 63 hinge, which has the effect of 
increasing the stiffness and natural frequency of tail-rotor flapping. 

The equation for tail-rotor pitch from reference 4 is given by: 

“ '^ct “ ®ot tan(63> ( 59 ) 


where 6^^: is the commanded tail-rotor collective and 63 is the lag offset angle. 

Notice that the tail-rotor pitch is coupled to the coning angle, aof Since equa- 
tions ( 44 ) and ( 59 ) are now coupled, they should be solved simultaneously. Eliminat 
ing between these two equations yields; 


0 


0 1 


®ct ■ *^3 

1 + Z2 tan 63 


( 60 ) 


where Zj^, Z^ , and Z3 are defined from equations (^ 4 ) by: 


)t ^i''c ^ ^ 3 ^lt 


( 61 ) 


Engine Model 

The SH- 3 G uses two T 38 -GE- 8 F gas turbine engines operating together except, of 
course, in the event of engine failure. Each engine has a gas-generator section pro- 
viding compressed air for a free or power— turbine section. An engine governor con- 
trols fuel flow to maintain constant power-turbine speed under changing loads. Fuel 
flow primarily affects the gas-generator speed which controls the torque applied to 
the power turbine. 

The governor is limited in a complex way during eri-»ine acceleration or decelera- 
tion to avoid the following undesirable conditions: (1) turbine overtemperature, 

( 2 ) compressor stall, ( 3 ) overrich flameout, or (A) overlean flameout. These limits 
are in direct opposition to obtaining maximum power and immediate response. The net 
result to the pilot is: ". . . engine response to new power or speed settings is not 
instantaneous: a few seconds must be allowed for the engine to stabilize at the new 

condition" (ref. 11, p. 1-10). 

A very important engine characteristic to simulate is this delay to sudden 
changes in the power required. A simple and effective way to model the engine is as 
a toniuc device regulating main rotor rpm. The main rotor is then modeled approxi- 
mately as a pure inertia and the controller is .in application of pseudo-derivative 
feedb.ick in a simple form (ref. 12). 

Figure 10 shows a block di.agr.am of the engine/governor model. As developed in 
refi'rence 12, the characteristic equation is: 
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Figure 10.- Engine and governor model. 


Is^ + K^s + K. = 0 (62) 

where d stands for derivative and i for integral. 

If equation (62) is critically damped, then the response to a step load» , is: 

This implies a characteristic time: 


1 



(64) 


Now from reference 8, the inertia of the main rotor and hub 10,190 slug-ft^ 
or I = 1067 f t“lb-sec/rpm. If t is chosen as 1.5 sec to simulate the engine delay, 
then = 1423 ft-*lb/rpm. Critical damping determines as 474 f t-lb/rpm-sec. 

The engine model is summarized by: 

ARPM = RPM ~ 203.3 \ 

J*RPM = J*RPM + ARPM x dT I ( 65 ) 

Qg = Ki J* RPM + K^(ARPM) ) 

where 


.\RPN deviation of main-rotor rpm from the nominal of 203.3 
J*RPM integral of rpm error 
DT cycletime 

engine torcjue 
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Control .Systi'.m and Rigging 


The flight-control -system model consists of rigging iniormntion, the 

automatic stabilization equipment (ARK) model, and the barometric altitude hold mode. 
Rigging is the mechanical gain, Including coupling, between the pilot's controls and 
tlie motion of the swashplatc, i.e., collective stick (in.) to collective pitch (eg). 
The automatic control system has four channels of stabilization: the coll(:ctive, the 

pitch, the roll and the yaw channels. The SH-3C is attitude-st.nbilized in pitch an 
roll. The yaw channel is a heavy yaw damper plus heading hold control by microswitches 
on the rudder pedals. The barometric altitude hold mode is controlled by a separate 
switch from the ARE and drives the collective channel. 

The rigging constants are taken for the most part from reference 13, the SH-3G 
maintenance manual. The ARE is taken chiefly from reference 8, the Navy trainer math 
model with some clarification from reference 11, the Navy flight manual. 

The SH-3G has four pilot controls: collective stick, lateral stick, longitudinal 

stick, and rudder pedals which control the main rotor collective pitch, the lateral 
flapping angle, the longitudinal flapping angle, and the tail rotor collective, 
respectively. 

Each pilot control has physical travel limits as detailed in table 1. Each con- 
trolled parameter, i.e., tail-rotor collective, also has physical limits as shown in 
table 2. 


TABLE 1.- PHYSICAL LIMITS OF THE PILOT'S CONTROLS (in.) 


Control 

Lower limit 

Collective, stick 

0.00 

Lateral cyclic stick 

-7.00 

Longitudinal cyclic stick 

-7.54 

Rotary rudder pedals 

-3.25 


Upper limit 


7.46 
7.00 

6 . 46 
3.25 


Sign convention 


+ UP 
+ RIGHT 
+ AFT 

+ NOSE RIGHT 


TABLE 2.- PHYSICAL LIMITS OF THE MAIN ROTOR AND TAIL ROTOR (deg) 


Rotor parameter 

Lower limit 

Upper limit 

Sign convention 

Main rtHor collective pitch 

8.10 

19.50 

+ UP 

Lateral flapping angle 

-9.10 

6.90 

+ RIGHT 

Longitudinal flapping angle 

-15.35 

10.15 

+ FWD 

Tail rotor collective pitch 

-6.50 

25.00 

+ NOSE LEFT 
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The form of the rlggioR oquntlona arc shown below; 


‘is - ‘^ln=‘u + '^U-e’^c 

®i8 ■ ■'la^la '^lo’'lo ‘^♦o ■" ''lo-e’‘c 


( 66 ) 


where Xc, Xia, Xjq, and Xj- are the collective stick, the lateral cyclic-stick, the 
longitudinal cyclic-stick, and the rudder-pedal positions, respectively, in inches. 
Note the coupling between collective and the other three controls. The cyclic control 
phase angle il/o, resulting from the rotor hinge offset, has been calculated by equa- 
tions in appendix H of reference 14. The rigging constants are shown in table 3. 

Each of the ASE components in equations (66) are limited to 10% control authority 
as shown in table 4: 


TABLE 3.- RIGGING 
CONSTANTS 


Constant 

Value 

Units 


2.72 

deg 

, 0 

8. 10 

deg 


9.25 

deg 

n 

C’” 

1.528 

deg/in. 

Kla 

1. 143 

deg/in. 

Klo 

-1.821 

deg/in. 

Kt 

-4.846 

deg/in. 

Kc"la 

-0. 1475 

deg/in. 

^C“lo 

-0.3485 

deg/in . 

K,-t 

1.135 

deg/in . 


TABLE 4.- ASE AUTHORITY 


Channel 

Limits, deg 

Symbol 

Collective 

±1.14 

ASE, 

Lateral 

±1.60 

^^^la 

Longitudinal 

±2.55 

ASE]_o 

Yaw 

±3.25 

ASE^ 


The control system equations are shown below: 

ASE, = tGj,^^(h - h^^j) + 

G, X, \ S. 


ASE 


la 


/ n A \ ^ase 


la" 

S 


( 67 ) 


ase 


- (fi„" + t'-pQ) r Ti 


ASE = 


G,.R + 0.(0 - 
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who.ro 

^bah barometric altitude hold switch (1 « on) 

^ase i^utomatic stabillMtlon switch (1 « on) 

'Vef barometric altitude reference: set when the barometric altitude hold switch is 

on 

^^la-’^la’^lo *^t time constants for first order filters 

’^ref heading reference: set when feet off the rudder pedals 

P*R body axis and yaw rates, respectively 

s Laplace Transform variable 

The control system gains and time constants are shown in table 5: 


r^NvVri 

c'f vl’’-' ** 


TABLE 5.- CONTROL SYSTEM GAINS 
AND TIME CONSTANTS 


Constant 

Value 

Units 

■^^la 

0.8 

sec 

"^la 

.625 

sec 

"^lo 

.625 

sec 

■^t 

.3125 

sec 

Gc 

.1430 

deg/in . 

Gbah 

-.01238 

deg/ft 

Go. 

-.1006 

deg/deg 

s 

-.1187 

deg/(deg/sec) 

^^la 

1.515 

deg/in. 

Go 

.2401 

deg/deg 

Gq 

.2593 

deg/ (deg/sec) 

^lo 

-.4504 

deg/in. 

Gr 

.2428 

deg/ (deg/sec) 

n 

H- 

1.3 

deg/deg 


The lateral and longitudinal channels are filtered to prevent sharp transients 
when the ASE is switched on. 


The actuator dynamics have been neglected in the control model because of the 
relatively short time constants involved (<0.2 sec). 


TRIM METHOD 


Ihe trim algorithm used is adapted from BQUIET, a general trimming subroutine for 
simulations at NASA Ames and documented in reference 15. BQUIET nulls six states with 
six or less contrtils by: (1) finding perturbations in the six states for each 
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control. (2) extrapolntlnp. linearly to find a new control vector which will minimize 
the. Mtate vector in the least-aquares sense, and (3) iterating until the. stares are 
within some tolerance of zero. Recognizing that the math model to be trimmed may be 
hlghi; mmlinear. (2) is modified by a ••gradient gain^^ of 0.5 i.e.. the extrapolation 
is only carried out half way. This damps oscillations about the solution and can 
prevent divergence . 

The BQUIET algorithm has been simplified and tailored for the SH-3G simulation. 
This required initializing the program properly between perturbations, choosing the 
controls and states, and making certain improvements in the details of the algorithm. 


Trim Initialization 


The SH-3G trim algorithm requires an estimate of the partial derivative of each 
state with respect to each control on every iteration. The partials are 
choosing a reference control vector, evaluating the corresponding reference state 
vector, disturbing each control in turn by a small percentage J-ts travel, an 
evaluating the chLge in the state vector from the reference. The state vector is 
evaluated by setting initial conditions of the math model (including the controls), 
allowing the model to "fly^' for two shortened cycles, and observing the states. The 
trlrc^le timfis chosen'very short so that integrals present in the math model will 
effectLely not operate in trim mode, avoiding the need for special loops around each 
integrator in triL Filters, however, must be set to their steady-state-gain value 
while trimming or they will be unaffected by changes in the controls. 

Ordinarily, the linear accelerations along the body axes 

erations about the body axes are chosen as states to be trimmed. All position, 

velocity, attitude, and angular velocity variables must therefore be set at the 
desired^trim values before each pass through the math model to determine the ^ccele - 
ations. Since, in this math model, two of the controls are attitude variables, this 
initialization must be recalculated before every pass. 

The initial velocities and attitudes- need to be specified in body axes well 
as in inertial axes. For the convenience of the researcher, sideslip angle, flight- 
path angle, the equivalent velocity (knots), and the wind vector .should be specifia- 
ble. The bank angle and angle of attack are appropriate control variables and are, 
therefore, also specified. The problem may be stated formally as: given, a, 6, Yv. 

- - ... r. ,, ,, .,. ggg appendix A for the solution. 


W, V„, 4>, find; V^, V, 0, i!<. 


Th> 

In addition, the model should be trimmable in hover, rearward flight, vertical 
climb and descent, and sideward flight. These problems are also adaressed in 
appendix A. 


Trim States and Controls 


The classical Newton’s method of finding the zeroes of a function consists 
finding the slope of the function at a trial point, extrapolating li^^ar^y to zero, 
evaluating the function at the new trial point, and repeating. For this application, 
this method is extended to many dimensions and the extent of extrapolation s con- 
trolled by a "gradient" gain. 


The original algorithm from BQUIET assumes an 
tlons, l.e., fewer controls than states. To avoid 


overdetermined set of linear equa- 
the possibility of uncontrollable 
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stntGfi and to simplify the computation, the present trim method insists on as many 
controls as states. This results in a square matrix of partial derivatives of the 
states with respect to the controls. Mathematically: 

aXi 

X - Ju where J, . = t — 

i.l 3u^ 

The SH-3G math model originally trimmed six states consisting of the i in ^ar 
accelerations along the body axes and the angular accelerations about the -body axes. 
The six controls were? the initial angle of attack, the iriitial bank angle, longi 
tudinal cyclic stick, lateral cyclic stick, rotary rudder-pedal position, and collec- 
tive stick. Trimming improved dramatically with the addition of the main-rotor angu- 
lar acceleration as a seventh state and the initial value of the engine model filter 
as a seventh control. 

A trim iteration, then, contains eight passes through the math model. Each of 
the first seven passes fills one column of the square seven by seven matrix of 
partial derivatives , the Jacobian . Each column of the Jacobian corresponds to per- 
turbations in the states due to a 0.01% of travel change in one control. After the 
seventh pass, the Jacobian is inverted and multiplied by the current reference state 
times the gradient gain to produce a new reference control vector. The eighth pass 
through the math model determines the new reference state vector which should be 
closer to trim. If the states are less than some trim criteria, the process is 
stopped . Mathematically : 

6u = -GJ‘“^X ; until |x^l < ; 0 < G < 1 


Trim Algorithm Improvements 

Three minor improvements have been made in the BQUIET algorithm « The improve- 
ments involve the limits on control changes, control restrictions when near the 
boundaries of their travel, and variations in the gradient gain. With these improve- 
ments the final math model trims at any airspeed from -30 knots to 135 knots and is 
fairly insensitive to the initial control guess. 

Some limit must be placed on the allowable change in each control, otherwise the 
trim process may converge very slowly or not at all. The SH-3G trim evaluates the 
new proposed change in the control vector, finds the control with the largest change 
as a percentage of its travel, limits this control to 10%, and re-scales the other 
controls to preserve the direction of the new control vector. 


If a control is allowed to exceed its travel, highly nonlinear response to per- 
turbations is likely because of physical limits and mechanical stops modeled. The 
SH-3G trim, after limiting the change in the control vector as above, makes a further 
check to see if any control would exceed its travel. If so, that control is limited 
to half the remaining travel and the control change vector is re-scaled. This is 
similar to what happens in American football when a 15 yard penalty is called against 
the offense and less than that remains to the goal line. This allows trimming arbi- 
trarily close to the control limits. 

The gradient gain, mentioned earlier, is a measure of the extent of extrapolation 
of the Jacobian matrix calculated on each iteration in trim. The ad hoc value for 
the gradient gain used in HQUTET is 0.5. Occasionally, this value is too large, 
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MODHl. VAhlOATtON 


The S11-3C math model validation Is ‘"s^t^^ariour^” 
lected to valldnte thoir own trainer simulation oi the ^ ^ . 
trlLmed control positions and time histories were selected to check the S11-3C math 

model fidelity both statically and dynamically. 


Static Cheeks 


lorq„o,. .ma tri-aa control po.,itlono ore vollclatcd for airspeeds vorylng 

fro,. -3 Unots o .U3^.nots.^^X^ 

tl.^is\ollLtod wro from 16.000- to 20,000-11. gross weigl.t nod sen level to 2,000 ft 
altitude. The results are presented in figures 11 to 13. 


Plpure 11 shows engiiu. torqne versus olrspeed for level ^ 

40 knots to 135 knots. A eo,.porlson I;.;-'™ the “ ' 
Shown for the medium gross weight ol 17. Tbs lb. ioi aiispecus incweci 



iv’uie 11.- level I 1 i gbt i>etf ormance . 



\?.{) knots t lu' torqiu's am almost liliMUloal, but tlu‘ math iiuulol is optimistically low 
below 70 knots ami above \2i) knots. At AO knots t lu» tlisi’ri'paney U\ larp.eJ^'l^ * proliably 
bissuise ol a mismat c'h in lilaile profile Orap,. 

Low speed and hover perfin'iiianee is slu^wn In fipuri' 12. Indicated torque for 
airspeeds from U) kiuUs aft or rearward Ui AO knots forward is plotted for bi>th tlie 
math model and t lu' fUpht ti'st. As can lu' seen, the flipht data are somewhat scat- 
tered. It doc’S appear, however, that the math model Is a few percent low especially 
in hovt'r and la.'arward f 1 I pht . 

I'Mpiires llUn) and I'Kb) compare the trimmed control positions for airs-peeds from 
30 knots to 130 knots in forward flipht. I’he sideslip anj7,le has been matched to the 
flight-test data for each airspeed. 

The collective position shows good agreement between flight test and math model. 
The minimum ct>lleetive position falls at 60 knots for the math model, which may be 
low depending on how the scatter in the flight data is read. This would explain, 
however, why the math model collective positions are high above 60 knots and low 
below 60 knots, i.e., a better match might be obtained by shifting the math model 
curve to the right, correspond! ng to an increase in induced power and a decrease in 
flat-plate drag power. 

'the lateral and longitudinal cyclic positions of the math model agree fairly well 
with the flight-test data, not varying more than about 7%. The longitudinal sensitiv- 
ity of the math model is slightly higher than the flight data. The math model lateral 
cyclic shows a decreasing trend wliile the flight data Increase, although the downward 
curvature of the two curves is similar. 
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Figure 13. - Concluded. 

The math model trimmed pedal position shows the correct trend but the wrong 
magnitude; being low by as much as 10%. This corresponds to a pedal displacement of 
5/8 of an inch, so it is questionable whether a pilot would notice the discrepancy. 

The pitch and roll attitudes show the correct trend. The trimmed roll attitudes 
are virtually indistinguishable between the math model and flight data; the pitch 
attitudes are a few degrees low. 

Figures 14(a) and (b) show the low speed and hover trimmed-control positions. 

Note that figures 13 and 14 are not comparable because of the large variation in gross 
weight. The effect of this can be seen in the collective position at 40 knots. In 
figures 13 the collective position is about 60% corresponding to the gross weight of 
19,017 lb as opposed to about 50% and 17,622 lb in figures 14. 

Again there is good agreement between the math model and flight data for the 
collective position and the roll and pitch attitudes. The lateral and longitudinal 
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Fij’uri* 14.- Low spood trimmod oontrol positions. 

I'vollo positions also show tho corrooL trends and do not difter from the fl 1 p.ht data 
hy more than alnnit 7".:. Tho longitudinal sensitivity of the math model is now less 
than t hi' tlijtht data in eontrast to lipaires 11. This is probably due to slip.ht dif- 
ferenees in venter of paavitv position whieh is fixed in the math model. 

I'lu’ math model trimmed pedal pi’sii ion is ap.ain low eompared to the t 1 i p,bf dat.i 
.ilthouy.h the trend appears eorreet . I’he pedal position is probably about 'ix low on 
averap.e as it was in tip.ures 11. This error is within reason, as the Sll-lt: ean be 
adjusted meehanieallv hv as mueh as 10... 
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Figure 13.- Static Interal-di rectii>nal stability. 

Overall, the SH-30 math model agrees very well with flight data in performance 
rimmed control positions throughout the airspeed range of the helicopter. 


Dynamic Checks 

ihe dynamic validation consists of a comparison with flight-test data for a 
step in aft cyclic, a 1 in. step in rlglit cyclic, a 1 in. step in right pedal, 
I ;10:. t.M-que increase tcollective step) .at .t trim airspeed of 70 knots. Nominal 
ht conditions were U»,00n-lh gross weight, ;',000 ft pressure altitude, and 
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Figure 15.- Concluded. 

standard temperature. In each case the pitch, roll, and yaw attitudes and pitch, 
roll, and yaw rates are examined. The results are presented in figures 16 to 19. 

Figure 16 compares flight test data for a 1 in. aft cyclic step with the math 
model response for the same 1 in. cyclic step. Both the pitch rate and pitch atti- 
tudes are well modeled with approximately the right shape and amplitude. The roll 
and yaw parameters show the correct trend toward the end of the run and the math 
model is lightly coupled as is the actual helicopter. 

A I in. right cyclic step is shown in I igure 17. As can be seen, the shape in 
roll and roll attitude are correct and the magnitude of the roll attitude is approxi 
mately correct. The roll rate magnitude, however, is low for the math model bv a 
factoi ol reaching 20 /sec in the flight data and only 10* /sec for the math model 
The effect of this can be seen in the roll response: the math model requires 3 

to reach maximum attitude-, while the flight data show only about 2 sec to reach 
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Figure 16.- Oae inch aft cyclic step. 

maximum attitude change. This discrepancy remains unexplained, as the math model 
will not roll at liO”/sec with a 1 in. right cyclic step even with tUe control system 
off and only the main rotor natural damping operating. 

The yaw attitude and yaw rate are well modeled showing the right shape and magni' 
tude. The pitch attitude and rate is lightly coupled in both the math model and the 
flight data . 

Figure 18 shows the response to a 1 in. right pedal stop. The yaw rate and yaw 
attitude are very well modeled in shape and magnitude. The roll attitude is also 
well modeled in sliape and magnitude. The math model roll rate is about double the 
flight-test data, reaching 7" /sec as opposed to 2‘’/sec in the flight data. The shape 
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Figure 17.- One inch right cyclic step. 

of the math modeled, roll-rate curve is correct, however. The pitch attitude and 
rate is very well modeled, although coupling is low. 

A 20% torque increase or, effectively, a collective step is shown in figure 19 
Ihe yaw and yaw rate are particularly well modeled with similar frequency and ampli 
tilde between the 1 light data and math model. The roll rate and roll attitudes of t 
math model seem to be at a lower frequency and amplitude than the flight data. The 
pitch and pitch rate show rouglily the same amplitude and shape for both the math 
model and flight-test data. 
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Figure 18.- One inch right pedal step. 


Overall the dynamic response in the primary axes, e.g., pitch response for 

ts^,.o.lcU..l vcr, ..n .xcent for tl.o roll rato. The coupled 

response Ls modeled fairly well and is small in any event. 


CONCl.USlONS 


A ,1 1 1 ,r th.. Sikorskv Sll-'Kl heliconter well-matched to flight-test data 

,„d l:ff-H..c and rcal-U»c »l»ulaUon ha.s been developed at A.es 

R, ■1-arch Center. The model eont.ilns e<luaU..mi .'f motion, an atmospheric m.idt 1 , . 


PA - 2000 ft 
OAT = 4“C 


— FLIGHT DATA 

— MATH MODEL 





TIME, sec 


Figure 19.- 20% torque increase, 

Dryden wind model, the fuselage aerodynamics, nonlinear main-rotor and tail-rotor 
models, a second-order engine model, and the SH-3G control system and rigging. An 
improved trim algorithm has ;ilso been developed. 

1. The SH-3n math model performance and t rimmed-control positions for airspeeds 
from -30 to 135 knots agree well in trend and magnitude with Navy flight-test data. 

2. The math model dynamic response in attitude and rate to a 1 in. step in 
longitudinal cyclic, lateral cyclic, or pedals agrees fairly well in the primary axis 
with Navy flight-test data. 


3 . T h 0 fa 1 lowing a cl hoc p r o c c d u r o h a a improve^ d t: i ^ v. 9 1 / 1 1 i c t r :1m a c c u r a c.y of a 

implified holicoptor math modol: 

a. Rnduce the blade lift curve alope until the tr;lmmod collective poaition 
in hover matchea flight data. 

b. Increase the induced power term in the main rotor torque equation until 
engine torque in hover matches flight data. 

c . Increase the fuselage f la t-plate drag area until engine torque at 
90 knots matches flight data. 
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APPENDIX A 


ORHSINAi. 

OF POOR 


VELOCITY INITIAIJZATION 


Throe velocity vectors (air» wind, and inertial velocity) need to be initialized 
in two coordinate systems — earth frame and body frame. This appendix contains a 
description of the various coordinate systems used, the problem, the solution, and 
some comments about special cases. 

Earth axes relate to a triad of orthonormal vectors; D, where N points 

north, E points east, and D points down into the earth. Body axes relate to 
another triad of orthonormal vectors k, 9, where k points out of the nose of an 
aircraft. ^ points out of the right wing, and k points down perpendicular to the 
plane of x and y. A useful set of coordinates for this initialization problem is 
the u, m, n set or "path" coordinates; u points along the earth relative velocity, 
m is perpendicular to u and points right in the horizontal plane, and n points 
down perpendicular to the plane of u and m. These coordinate systems are summarized 
by the three, equations: 

V = V„N + V„E + V_D "earth" coordinates 

e N ED 

V =Vx+Vy+V£ "body" coordinates 

b X y-^ z 


V =Vu+Vm+Vn "path" coordinates 
p u m n 

The trim initialization problem can be formally stated as follows; 
Given : 

V the air velocity magnitude 

a 

wind component north 

W„ wind component east 

£ 

wind component down 

a angle of attack of the fuselage 

(S sideslip angle of the fuselage 

vertical flightpath angle 

Y, horizontal flightpath angle 

h 

({) Euler roll angle 


Find : 

u X body axis component of the air velocity vector 
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w. 




D 

Va 

Vaj 

Vai 

0 


N 


y body^fixlfi component: of the air velocity vector 
z body axis component of the air velocity vector 
north component of the earth relative velocity vector 
east component of the earth relative velocity vector 

down component of the earth relative velocity vector 

north component of the air velocity vector 
east component of the air velocity vector 

down component of the air velocity vector 

Euler pitch angle 


Euler heading angle 
Read the subscripts as follows: 

air h horizontal 

body N north 

wind E east 

v vertical D down 

The solution proceeds as follows: 

From the definition of angle-of-attack and sideslip calculate the body axis 
components of airspeed. 

u. = V C3Ca 
b a 

V. V Sb 

b a 

W, = V CbScx 

b a 

From the definition of airspeed, wind, and ground speed: 

V = V + W or 
a 

VG = V V + Ww 
a a 


(Al) 


Dot (Al) with u, m, and ri. 


V = V (V • u) + W(w • u) 
a a 


v„ + w 

u 


0 = V^(V.^ • m) + W(w • m) = 


m 


m 


0 = V (V • n) 4- W(w • n) = V. + W 


a a 


(A2a) 

(A2b) 

(A2c) 
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Hy definition: 


a 

■''au 

W'‘ 

- W' ■ 
n 
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'am "a. 




n 


(A 3b) 


Rewriting (A2b) and (A2c) 




= -W 


m 


-W 


(A2b) 

(A2c) 


Substituting (A2b) and (A2c) into (A3a): 


= Va + W" + 

a ^u m n 


Solving (A3b) for (W^ + W„) and substituting into (A4): 


Substituting into (A2a): 


vi = - Wf ) 

a u 


V » /V? - (W^ - wO + w 

a u u 


(A4) 


(A5) 


Transforming a vector from earth coordinates to path coordinates requires two pure 
rotations . 

W = Ty, W 

p ' V ' h e 

where 


Vv 



0 

-Sy 

’ V 




0 

0 

1 

0 



^^h 

0 

Sy 

’ V 

0 

C Y 

V 


0 

0 

1 


Applying the above transformation: 



b 



And the component of tlie wind parallel to the earth relative velocity vector is: 


W, " '''v'-'h • ''N ^ ' 


W.. - S-, 


w. 


(A6) 
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OF POOR QUALITY 

With the earth relative velocity magnitude from equation (A5) and the definition 
Yv «nd Y)i: 

Ve = V Cy^SY^ 

Vd = V(-Sy^) 

Equation (Al) now gives: 



From the definition of body axes: 



where 0 and are still unknown and T^, Tq, and are defined in equation ( 
the main text. 

Separating 9 and ip: 



Note that (A7b) is independent of 0 and (A7c) is independent of t|/: 

C<t) - S(^ = + ^'a£ C 4 ) 

-U^ S'; + (V^ ScJ) + C<1>)C0 = 
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(A7a) 

(A7b) 

(A7c) 

(A7b) 

(A7c) 


OWG<NAU PAGE W 
OF POOR 

Equations (A7b) and (A7c) can now be put in the form: 

A ain(q) B cos(q) 4- C 
This can be solved for q as follovrs: 


sin(q) 

Ls/K ^ f 




cos(q) 


cos(r)sin(q) - sin(r)cos(q) = 


where tan(r) = B/A; then sin(q-r) = 
for r : 


T* Sol\'ing for q and substituting 


q = arctan(~j + arcsinf 


Rewriting (A7b) in the form of (A8) : 


S^p = Vag C.(' + (W^ S<t. - Zi,) 


Applying (A9) : 


!;■ = arctan 


/Vap\ S+ - Vp 


where 


'ag = Pk + * 0 and Va^ > 0 


Rewriting (A7c) in the form of (A8): 


Uu so = (V, S 4 ^ 4 w, C()))cn - 


Applying (A9) : 


arctan I — ,1 I + arcsin 

V ' «1ji' + (V. S(J) + W, C^)'2 

^ U D D 


where lb ^ 0 and V * 0. 
b a 

Note that for c*) = i = 0: 


. - arctan 


m = . 
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as required. Rearward flight is handled by adding 180° to 3 and Yy. This— implies 
(X and should be limited to ±90°. This limitation results from the 
definitions : 


V = Ty Tv, V 
p 'V 'h e 

In words, the air-velocity direction is defined by two rotations (a and 3) from 
body axes and the path-velocity vector is defined by two rotations (y^^ and Yv) fi’oai 
earth axes. To rotate the velocity vector to the opposite direction, the last rota- 
tion angle (3 or Yv) ®ust be used, instead of a or Yh» the resulting vector will 
not, in general, point in the opposite direction. 

Hover is handled by not allowing the total air velocity to equal zero exactly. 

If the air-velocity magnitude is less than some small criterion, then the air velocity 
is assigned the criterion value times the sign of the velocity. Mathematically: 

If l\i “ sign(Vg) x 

The method described in this appendix is valuable for solving other similar 
problems. Particularly useful are the two "tricks" embodied in equations (A5) 
and (A9). 
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appendix b 


SYMBOLS 


ASE^ 

ASE, 

la 

ASEt 

lo 

ASEj. 

a 

a' 


“ot 

ai.b^ 


3 1 , b 1 


B 

b 


lateral and longitudinal cyclic control angles in shaft axes 
collective channe.' automatic-control command 
lateral channel automatic-control command 
longitudinal channel automatic-control command 

yaw channel automatic-control command 

blade lift curve slope, rad 

angle of tilting of the thrust vector due to lift 
rotor coning angle 

tail rotor blade coning angle 

lateral and longitudinal flapping angles in control axes 
lateral and longitudinal flapping angles in shaft axes 
blade tip loss factor 
number of blades 


• • • ^10 moment of inertia coefficients (eq. (12)) 
fuselate flat-plate drag area 

CFd, fuselage drag-force coefficient due to local angle of attack 

^'^^2 fuselage drag-force coefficient due to local sideslip 

CFli fuselage lift-force coefficient due to local angle of attack 

fuselage lift-force coefficient due to sideslip 
^Fy fuselage side-force coefficient 

foil ing-moment coefficient 

pitching-moment coefficient due to local angle of attack 
pitching-moment coefficient due to sideslip 
yawing-moment coefficient 
lorciue coefficient, Q/i- " :R'’ 
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ct, 


m 


D 

DT 

e 

^kf 


F Q 9 ^ A * ^ 3 . 


F^’^y >^z 

^Xs»Pys»^^Zj 

G 

G„ 


bah 


Gx 


^la 

lo 


G,, 

c; 


ref 


page is 

qUAUTV t .. 

main rot: or- thrust coefficient 

rotor side-force coefficient 

rotor blade chord 

fuselage drag force, (+ aft) 

computer cycle time or integration step 

rotor flapping-hinge offset 

fuselage angle-of-attack correction for main rotor downwash 

main rotor downwash factor 

fuselage aerodynamic force components 

gravitational force 

force in earth axes 

force in body axes 

force in sh-ift axes 

gradient gain on the Jacobian matrix 
control gain on yaw rate 
control gain on altitude error 
control gain on collective stick position 
control gain on roll rate 

gear ratio between the main and tail rotors, 6*28937 

control gain on lateral cyclic-stick position 

control gain on longitudinal cyclic-stick position 

control gain on pitch angle 

control gain on bank angle 

control gain on heading error 

rotor drag force 

altitude above sea level 

control reference altitude, set when the altitude hold switch is 
turned on 


OWGtNAU 
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I 

T T T 

J 

K. 


K 


c-la 


K 


c-lo 


K 


c-t 


K. 


K, 


K 


la 


K 


lo 


'sw 


K 


nu 

Kt 

L.M.N 

Lo 

LAT 

LON 




M 

m 

"Htf 

n 

I’,Q,R 

P,T 


effective inertia of the main rotor, -Ijn 

moment of inertia of the main rotor blades and hub 
# 

body axis moments of inertia 
rotor side force 

rigging constant* between collective stick and collective pitch 

coupling between collective stick and lateral swashplate angle 

coupling between collective stick and longitudinal swashplate angle 

coupling between collective stick and^.tail-rotor pitch 

engine governor constant affecting the derivative of Arpm in the 
characteristic equation 

engine governor gain on the time integral of rpm error 

rigging constant between lateral cyclic stick and lateral cyclic-control 
angle 

rigging constant between longitudinal cyclic stick and longitudinal 
cyclic-control angle 

thrust correction for rotor spanwise flow not accurately modeled 
thrust correction for nonuniform flow 

rigging constant between pedal position and tail-rotor collective 

rolling, pitching, and yawing moments in body axes 

magnitude of a generalized step load on the engine 

latitude of the helicopter 

longitude of the helicopter 

fuselage lift force (+ up) 

pure rolling, pitching, and yawing moments in shaft axes 
Mach number 

mass of the helicopter 

mass moment of a rotor blade 

count of current time step 

body axes roll, pitch, and yaw rates 

free-stream temperature and pressure 


A 7 


0 
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. 'It ✓ 


*^c 

Pg.Qg.Rjj 

Qe 

Qm 




ase 


bah 


s 

T 


Ta£»Ta„»Ta^ 


"m “n 


Ta,,Pa, 


Tr.Tr 

Uc.Vc.Wc 

Ug.Vg.Wg 

V 

Va 


fir 


VaN’VaE»V,,D 


Va ,V.., ,Va^ 

u m n 


cn 1 


' eq 


control axes roll, pitch, and yaw rates 

shaft axes roll, pitch, and yaw rates 

engine torque 

main-rotor torque 

tail-rotor torque 

dynamic pressure 

compressible dynamic pressure, as measured by a pitot static tube 
rotor disc radius 

radius of the earth, 20,898,908 ft 

automatic stabilization equipment switch (1 = On) 

barometric altitude hold switch (1 = On) 

Laplace transform variable 
rotor thrust 

fuselage aerodynamic moments in body axes 

ratio of standard temperature and pressure at altitude to sea-level 
values 

total temperature and pressure ratios 

total temperature and pressure 

airspeed in body axes 

airspeed in control axes 

airspeed in shaft axes 

earth relative velocity magnitude 

airspeed magnitude 

ground speed 

airspeed in earth axes 

airspeed in path axes 

calibrated airspeed corrected for the effects of compressibility 
equivalent airspeed at sea- level density 
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X, 


c 

Xxa 

X, 


0 


‘^2 >^3 

a 


3 

3r 

Y 

■^h 

Y, 


V 

ARPM 


AT 

6 

^0 > 1 ’ ^ 2 

6, 


Or 


'ot 


Ot , 0 




earth relative velocity In earth axes 
wind velocity In earth axea 
wind velocity In path axes 
collective stick position 

coordinates of the rotor hub relative to the CG in body axes 

lateral cyclic^stick position 

longitudinal cyclic-stick position 

tail rotor pedal position 

fuselage side force (+ right) 

rotor constants defined by equation (61) 

fuselage centerline angle of attack 

fuselage centerline angle of attack corrected for main rotor downwash 

average rotor-blade angle of attack in hover 

sideslip angle of fuselage (+ nose left of velocity vector) 

rotor orientation angle 

rotor lock number 

horizontal flightpath angle (eq. (5)) 
vertical flightpath angle (eq. (5)) 
main-rotor rpm deviation from nominal (203.3) 
tempera ture above standard conditions 
rotor-blade profile-drag coefficient 

coeff .ients of the Bailey profile-drag coefficient (eq. (49)) 

tail-rotor 6^ hinge phase angle 

rotor collective pitch angle 

tail rotor collective pitch angle 

lowest main-rotor collective pitch 

neutral tail-rotor collective pitch with collective in full low position 
rotor blade twist from root to tip 
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'7C 


«ct 


'tn 




la 




la 




(f)g , 0g 
^0 




ref 


Vectors 

N,E,D 

u,ni,n 

5^ $ y j 2 

x,u 

t 

Matrices 


T.,T 


blade pitch at 75% of the rotor-disc radius 
tail-rotor commanded collective- pitch 


<r 


main-rotor inflow ratio 

main-rotor tip-speed ratio 

induced inflow ratio 

atmospheric density 

rotor solidity ratio 

lateral channel filter-time constant 

yaw channel filter- time constant 

lateral stick filter-time constant 

longitudinal stick filter-time constant 

lag time constant for changes in the inflow ratio 

roll, pitch, and yaw Euler angles 

rotor-shaft roll and pitch relative to the (-Z) body axis 

rotor flapping phase angle 

reference heading 

main rotor angular rate 

angular rate of the earth, 0.000072722 rad/sec 


unit vectors along earth axes 

unit vectors along path axes 

unit vectors along body axes 

state and control vectors for trimming 

trim-error criteria for stopping the trim process 


Jacobian matrix-matrix of partial derivatives 

rotation matrices for locating the airspeed vector relative to 

axes 




n- ,Xy, rotation matrices for locating the earth relative vector on the earth 

'V 'h 

T ,T, ,T rotation matrices for converting earth axes to body axes 

ij) 0 i|i 

Sub scripts 

b body axes 

c control axes (aligned to the axis of no feathering) 

e earth axes 

m main rotor parameters 

p path coordinates (see appendix A) 

s shaft axes (aligned to the rotor shaft axis) 

t tail-rotor parameters 

(’) 


time derivative 
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